Introduction

< This is a preliminary OM that should be subject to revision >

This operating model was specified using a stock assessment document (referred to as the “Assessment” herein):

DFO. 2012. Stock Assessment for the inside population of Yelloweye Rockfish, (Sebastes ruberrimus), In British Columbia, Canada for 2010. DFO Can. Sci. Advis. Sec. Sci. Advis. Rep. 2011/084 13p

and a COSEWIC report (referred to as the “COSEWIC report” herein):

COSEWIC. 2008. COSEWIC assessment and status report on the Yelloweye Rockfish Sebastes ruberrimus, Pacific Ocean inside waters population and Pacific Ocean outside waters population, in Canada. Committee on the Status of Endangered Wildlife in Canada. Ottawa. vii + 75 pp.(www.sararegistry.gc.ca/status/status_e.cfm)

Following the provision of catch composition data, mean length data and indices, the operating model should be conditioned using Stochastic SRA. This will provide estimates of recruitment deviations, stock depletion and selectivity.

Map courtesy of Sean Anderson’s groundfish synopsis

Robustness tests

No robustness OMs have been produced yet but could address uncertainties relating to:

  • the relative abundance index used to condition the OM

  • the selectivity of fishing

  • previous time-area closures

Operating Model

The OM rdata file can be downloaded from here

Download and import into R using myOM <- readRDS('OM.rdata')

Species Information

Species: Sebastes ruberrimus

Common Name: Yelloweye Rockfish

Management Agency: DFO

Region: British Columbia

Latitude: 54.5, 53, 51.2, 48.2, 48.2, 50, 54.5

Longitude: -134, -133.5, -131, -125.8, -122, -122, 129

OM Parameters

OM Name: Name of the operating model: Canary_Rockfish_BC_DFO

nsim: The number of simulations: 192

proyears: The number of projected years: 50

interval: The assessment interval - how often would you like to update the management system? 4

pstar: The percentile of the sample of the management recommendation for each method: 0.5

maxF: Maximum instantaneous fishing mortality rate that may be simulated for any given age class: 0.8

reps: Number of samples of the management recommendation for each method. Note that when this is set to 1, the mean value of the data inputs is used. 1

Source: A reference to a website or article from which parameters were taken to define the operating model

Stock assessment for Canary rockfish (Sebastes pinniger) in British Columbia waters. Stanley R. Starr P. Olsen N.

Custom Parameters

Stock depletion (SSB(2009) / SSB(1918)) was prescribed as a log-normal random variable with mean 0.123 and CV 0.43.

Stock Parameters

Mortality and age: maxage, R0, M, Msd

maxage: The maximum age of individuals that is simulated. There are maxage+1 (recruitment to age-0) age classes in the storage matrices. maxage is the plus group where all age-classes > maxage are grouped, unless option switched off with $plusgroup=0 . Single value. Positive integer.

Specified Value(s): 110

The Assessment cites a maximum age of 110.

R0: Initial number of unfished recruits to age-0. This number is used to scale the size of the population to match catch or data, but does not affect any of the population dynamics unless the OM has been conditioned with data. As a result, for a data-limited fishery any number can be used for R0 . In data-rich stocks R0 may be estimated as part of a stock assessment, but for data limited stocks users can choose either an arbitrary number (say, 1000) or choose a number that produces simulated catches in recent historical years that are similar to real world catch data. Single value. Positive real number.

Specified Value(s): 1e+05

Unless management options are specified in absolute numbers (e.g. tonnes) the MSE is scale-less (has no units) and this value is inconsequential. Here it is set to 100,000 arbitrarily.

M: The instantaneous rate of natural mortality. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. Uniform distribution lower and upper bounds. Non-negative real numbers.

Specified in cpars: 0.02, 0.03

The COSEWIC report cites values in the range of 0.02-0.03

Msd: Inter-annual variation in M expressed as a coefficient of variation of a log-normal distribution. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. If this parameter is positive, yearly M is drawn from a log-normal distribution with a mean specified by log(M) drawn for that simulation and a standard deviation in log space specified by the value of Msd drawn for that simulation. Uniform distribution lower and upper bounds. Non-negative real numbers

Specified Value(s): 0.05, 0.1

I set an arbitrary level of inter-annual variability with a lognormal CV of between 5% and 10% (i.e. 0.05 to 0.1), corresponding with 95% probability interval of approximately +/-10% to +/- 20%. Note that due to the longevity of Yelloweye Rockfish, quite substantial inter-annual variability in M would be necessary to generate data inconsistent with the assumption of time-invariant M (noting the possible exception of a trend in M, below).

Natural Mortality Parameters

Sampled Parameters

Histograms of simulations of M, and Msd parameters, with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

The average natural mortality rate by year for adult fish for 3 simulations. The vertical dashed line indicates the end of the historical period:

M-at-Age

Natural mortality-at-age for 3 simulations in the first historical year, the last historical year (i.e., current year), and the last projected year:

M-at-Length

Natural mortality-at-length for 3 simulations in the first historical year, the last historical year (i.e., current year), and the last projected year:

Recruitment: h, SRrel, Perr, AC

h: Steepness of the stock recruit relationship. Steepness governs the proportion of unfished recruits produced when the stock is at 20% of the unfished population size. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This value is the same in all years of a given simulation. Uniform distribution lower and upper bounds. Values from 1/5 to 1.

Specified Value(s): 0.55, 0.7

Without stock-recruitment model estimates (or spawning data and a recruitment index), I borrow from the Canary Rockfish Assessment

SRrel: Type of stock-recruit relationship. Use 1 to select a Beverton Holt relationship, 2 to select a Ricker relationship. Single value. Integer

Specified Value(s): 1

A value of 1 represents the Beverton-Holt stock recruitment curve which is commonly assumed form of density dependence for rockfishes in BC.

Perr: Recruitment process error, which is defined as the standard deviation of the recruitment deviations in log space. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. Uniform distribution lower and upper bounds. Non-negative real numbers.

Specified Value(s): 0.25, 0.32

Estimated by Stochastic SRA (Walters et al. 2003)

AC: Autocorrelation in the recruitment deviations in log space. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided, and used to add lag-1 auto-correlation to the log recruitment deviations. Uniform distribution lower and upper bounds. Non-negative real numbers.

Specified Value(s): 0.75, 0.81

Estimated by Stochastic SRA

Recruitment Parameters

Sampled Parameters

Histograms of 48 simulations of steepness (h), recruitment process error (Perr) and auto-correlation (AC) for the Beverton-Holt stock-recruitment relationship, with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Time-series plot showing 3 samples of recruitment deviations for historical and projection years:

Growth: Linf, K, t0, LenCV, Ksd, Linfsd

Linf: The von Bertalanffy growth parameter Linf, which specifies the average maximum size that would reached by adult fish if they lived indefinitely. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This value is the same in all years unless Linfsd is a positive number. Uniform distribution lower and upper bounds. Positive real numbers.

Specified in cpars: 74.095, 78.985

The value for females is taken from Figure 13 of the COSEWIC report. Arbitrarily, a small degree of uncertainty is used to bracket the mean value of 76.41cm

K: The von Bertalanffy growth parameter k, which specifies the average rate of growth. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This value is the same in all years unless Ksd is a positive number. Uniform distribution lower and upper bounds. Positive real numbers.

Specified in cpars: 0.017, 0.081

Similarly to Linf the von B. k value of 0.02 (of the COSEWIC report) is bracketted by +/- 10%.

t0: The von Bertalanffy growth parameter t0, which specifies the theoretical age at a size 0. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. Uniform distribution lower and upper bounds. Non-positive real numbers.

Specified Value(s): -20.73, -20.73

Theoretical age at length zero. The value for females is taken exactly from the COSEWIC report and no uncertainty is simulated.

LenCV: The coefficient of variation (defined as the standard deviation divided by mean) of the length-at-age. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided to specify the distribution of observed length-at-age, and the CV of this distribution is constant for all age classes (i.e, standard deviation increases proportionally with the mean). Uniform distribution lower and upper bounds. Positive real numbers.

Specified Value(s): 0.1, 0.15

Arbitrarily assigned a CV in length at age of between 10 and 15 per cent following Figure 13 of the COSEWIC report.

Ksd: Inter-annual variation in K. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. If this parameter has a positive value, yearly K is drawn from a log-normal distribution with a mean specified by the value of K drawn for that simulation and a standard deviation (in log space) specified by the value of Ksd drawn for that simulation. Uniform distribution lower and upper bounds. Non-negative real numbers.

Specified Value(s): 0.2, 0.3

Studies of temporal variability in growth of rockfish are not common but older studies of rockfish growth have found moderate inter-annual variability in K values among years with a CV of around 25% (e.g. widow rockfish, Pearson and Hightower 1991).

Linfsd: Inter-annual variation in Linf. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. If this parameter has a positive value, yearly Linf is drawn from a log-normal distribution with a mean specified by the value of Linf drawn for that simulation and a standard deviation (in log space) specified by the value of Linfsd drawn for that simulation. Uniform distribution lower and upper bounds. Non-negative real numbers.

Specified Value(s): 0.05, 0.1

Taken from the same study above, this is relatively constant.

Growth Parameters

Sampled Parameters

Histograms of simulations of von Bertalanffy growth parameters Linf, K, and t0, and inter-annual variability in Linf and K (Linfsd and Ksd), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

The Linf and K parameters in each year for 3 simulations. The vertical dashed line indicates the end of the historical period:

Growth Curves

Sampled length-at-age curves for 3 simulations in the first historical year, the last historical year, and the last projection year.

Maturity: L50, L50_95

L50: Length at 50% maturity. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. The L50 and L50_95 parameters are converted to ages using the growth parameters provided and used to construct a logistic curve to determine the proportion of the population that is mature in each age class. Uniform distribution lower and upper bounds. Positive real numbers.

Specified in cpars: 33.095, 72.274

The COSEWIC report specifies 50% maturity at ages 16-20. Using the lower bound on K and Linf, the growth curve predicts a length of 35.7cm= 74(1-exp(-0.018 x (16+20.73))). Using the upper bound on K and Linf, the growth curve predicts a length of 28.1 cm = 79(1-exp(-0.022 x (20+20.73))). This is quite a naïve basis for guessing length at maturity since it does not account for aging error (compresses uncertainty) but then exaggerates uncertainty by using lower and upper bound pairs of the K and Linf growth parameters that are typically negatively correlated.

L50_95: Difference in lengths between 50% and 95% maturity. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. The value drawn is then added to the length at 50% maturity to determine the length at 95% maturity. This parameterization is used instead of specifying the size at 95 percent maturity to avoid situations where the value drawn for the size at 95% maturity is smaller than that at 50% maturity. The L50 and L50_95 parameters are converted to ages using the growth parameters provided and used to construct a logistic curve to determine the proportion of the population that is mature in each age class. Uniform distribution lower and upper bounds. Positive real numbers.

Specified Value(s): 7.16, 9.35

Following maturity curves of other rockfish species such as Canary Rockfish we assume that the length increment to 95% maturity is rougly 20% of the length at 50% maturity.

Maturity Parameters

Sampled Parameters

Histograms of simulations of L50 (length at 50% maturity), and L95 (length at 95% maturity), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Maturity at Age and Length

Maturity-at-age and -length for 3 simulations in the first historical year, the last historical year (i.e., current year), and the last projected year:

Stock depletion and Discard Mortality: D, Fdisc

D: Estimated current level of stock depletion, which is defined as the current spawning stock biomass divided by the unfished spawning stock biomass. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This parameter is used during model initialization to select a series of yearly historical recruitment values and fishing mortality rates that, based on the information provided, could have resulted in the specified depletion level in the simulated last historical year. Uniform distribution lower and upper bounds. Positive real numbers (typically < 1)

Specified in cpars: 0.034, 0.342

Table 1 from the Assessment provides a reference case stock depletion that we prescribe in custom parameters (med = 123, CV=0.43)

Fdisc: The instantaneous discard mortality rate the stock experiences when fished using the gear type specified in the corresponding fleet object and discarded. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. Uniform distribution lower and upper bounds. Non-negative real numbers.

Specified Value(s): 0.8, 1

Discard mortality rate was set to a conservative range of 80 per cent to 100 percent, consistent with the values published by the Pacific Management Council

Depletion and Discard Mortality

Sampled Parameters

Histograms of simulations of depletion (spawning biomass in the last historical year over average unfished spawning biomass; D) and the fraction of discarded fish that are killed by fishing mortality (Fdisc), with vertical colored lines indicating 3 randomly drawn values.

Length-weight conversion parameters: a, b

a: The alpha parameter in allometric length-weight relationship. Single value. Weight parameters are used to determine catch-at-age and population-at-age from the number of individuals in each age class and the length of each individual, which is drawn from a normal distribution determined by the Linf , K , t0 , and LenCV parameters. As a result, they function as a way to scale between numbers at age and biomass, and are not stochastic parameters. Single value. Positive real number.

Specified Value(s): 0

From the COSEWIC report Figure 12.

b: The beta parameter in allometric length-weight relationship. Single value. Weight parameters are used to determine catch-at-age and population-at-age from the number of individuals in each age class and the length of each individual, which is drawn from a normal distribution determine by the Linf , K , t0 , and LenCV parameters. As a result, they function as a way to scale between numbers at age and biomass, and are not stochastic parameters. Single value. Positive real number.

Specified Value(s): 3.32

From the COSEWIC report Figure 12.

Spatial distribution and movement: Size_area_1, Frac_area_1, Prob_staying

Size_area_1: The size of area 1 relative to area 2. The fraction of the unfished biomass in area 1. Please specify numbers between 0 and 1. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. For example, if Size_area_1 is 0.2, then 20% of the total area is allocated to area 1. Fishing can occur in both areas, or can be turned off in one area to simulate the effects of a no take marine reserve. Uniform distribution lower and upper bounds. Positive real numbers.

Specified Value(s): 0.1, 0.1

Population is uniformly distributed: size area 1 is same as frac area 1.

Frac_area_1: The fraction of the unfished biomass in area 1. Please specify numbers between 0 and 1. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. For example, if Frac_area_1 is 0.5, then 50% of the unfished biomass is allocated to area 1, regardless of the size of area 1 (i.e, size and fraction in each area determine the density of fish, which may impact fishing spatial targeting). In each time step recruits are allocated to each area based on the proportion specified in Frac_area_1. Uniform distribution lower and upper bounds. Positive real numbers.

Specified Value(s): 0.1, 0.1

Frac_area_1 is the default area for testing marine reserves or simulating habitat that is outside of the range of fishing. To simulate a mixed stock that may be subjec to a 10% MPA we simulate a stock in which 10% of individuals are in area 1 and 90% are in area 2, 0.1.

Prob_staying: The probability of individuals in area 1 remaining in area 1 over the course of one year. Please specify numbers between 0 and 1. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. For example, in an area with a Prob_staying value of 0.95 each fish has a 95% probability of staying in that area in each time step, and a 5% probability of moving to the other area. Uniform distribution lower and upper bounds. Positive fraction.

Specified Value(s): 0.85, 0.95

To simulate uncertain mixing among areas I assume that between 85 per cent and 95 per cent of individuals remain in the same area among years.

Spatial & Movement

Sampled Parameters

Histograms of 48 simulations of size of area 1 (Size_area_1), fraction of unfished biomass in area 1 (Frac_area_1), and the probability of staying in area 1 in a year (Frac_area_1), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Fleet Parameters

Historical years of fishing, spatial targeting: nyears, Spat_targ

nyears: The number of years for the historical simulation. Single value. For example, if the simulated population is assumed to be unfished in 1975 and this is the year you want to start your historical simulations, and the most recent year for which there is data available is 2019, then nyears equals 45.

Specified Value(s): 90

The data from the assessment run from 1920-2009, a total of 90 years.

Spat_targ: Distribution of fishing in relation to vulnerable biomass (VB) across areas. The distribution of fishing effort is proportional to VB^Spat_targ. Upper and lower bounds of a uniform distribution. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This parameter allows the user to model either avoidance or spatial targeting behavior by the fleet. If the parameter value is 1, fishing effort is allocated across areas in proportion to the population density of that area. Values below 1 simulate an avoidance behavior and values above 1 simulate a targeting behavior.

Specified Value(s): 1, 1

Spatial targetting. We’re going to stick to the default level of 1 (effort distributed in proportion to density)

Trend in historical fishing effort (exploitation rate), interannual variability in fishing effort: EffYears, EffLower, EffUpper, Esd

EffYears: Vector indicating the historical years where there is information available to infer the relative fishing effort expended.This vector is specified in terms of the position of the year in the vector rather than the calendar year. For example, say our simulation starts with an unfished stock in 1975,and the current year (the last year for which there is data available) is 2019. Then there are 45 historical years simulated, and EffYears should include numbers between 1 and 45. Note that there may not be information available for every historical year, especially for data poor fisheries. In these situations, the EffYears vector should include only the positions of the years for which there is information, and the vector may be shorter than the total number of simulated historical years (nyears).

The recent assessment predicts biomass from 1920 to 2009. Catches reported from 1950-2004 were used to calculate the implied F trend. These calculation can be found in the ‘Fcalcs’ worksheet of the Yelloweye_Rockfish_BC_DFO.xlsx workbook.

EffLower: Lower bound on relative fishing effort corresponding to EffYears. EffLower must be a vector that is the same length as EffYears describing how fishing effort has changed over time. Information on relative fishing effort can be entered in any units provided they are consistent across the entire vector because the data provided will be scaled to 1 (divided by the maximum number provided).

As above

EffUpper: Upper bound on relative fishing effort corresponding to EffYears. EffUpper must be a vector that is the same length as EffYears describing how fishing effort has changed over time. Information on relative fishing effort can be entered in any units provided they are consistent across the entire vector because the data provided will be scaled to 1 (divided by the maximum number provided).

As above

EffYears EffLower EffUpper
1920 0.00373 0.00373
1921 0.00380 0.00380
1922 0.00383 0.00383
1923 0.00383 0.00383
1924 0.00387 0.00387
1925 0.00385 0.00385
1926 0.00383 0.00383
1927 0.00387 0.00387
1928 0.00389 0.00389
1929 0.00391 0.00391
1930 0.00399 0.00399
1931 0.00399 0.00399
1932 0.00399 0.00399
1933 0.00403 0.00403
1934 0.00399 0.00399
1935 0.00397 0.00397
1936 0.00399 0.00399
1937 0.00399 0.00399
1938 0.00405 0.00405
1939 0.00413 0.00413
1940 0.00411 0.00411
1941 0.00409 0.00409
1942 0.00407 0.00407
1943 0.00411 0.00411
1944 0.00429 0.00429
1945 0.00453 0.00453
1946 0.00453 0.00453
1947 0.00453 0.00453
1948 0.00464 0.00464
1949 0.00472 0.00472
1950 0.00464 0.00464
1951 0.00939 0.00939
1952 0.00279 0.00279
1953 0.00593 0.00593
1954 0.00217 0.00217
1955 0.00351 0.00351
1956 0.00242 0.00242
1957 0.00223 0.00223
1958 0.00194 0.00194
1959 0.00469 0.00469
1960 0.00368 0.00368
1961 0.00337 0.00337
1962 0.00233 0.00233
1963 0.00375 0.00375
1964 0.00376 0.00376
1965 0.00429 0.00429
1966 0.00200 0.00200
1967 0.00337 0.00337
1968 0.00301 0.00301
1969 0.00308 0.00308
1970 0.00596 0.00596
1971 0.00590 0.00590
1972 0.00526 0.00526
1973 0.01540 0.01540
1974 0.00298 0.00298
1975 0.00240 0.00240
1976 0.00189 0.00189
1977 0.01090 0.01090
1978 0.01150 0.01150
1979 0.02920 0.02920
1980 0.01790 0.01790
1981 0.01390 0.01390
1982 0.00916 0.00916
1983 0.00731 0.00731
1984 0.01280 0.01280
1985 0.02420 0.02420
1986 0.03690 0.03690
1987 0.04870 0.04870
1988 0.06370 0.06370
1989 0.06640 0.06640
1990 0.07430 0.07430
1991 0.07050 0.07050
1992 0.02610 0.02610
1993 0.03780 0.03780
1994 0.08600 0.08600
1995 0.04050 0.04050
1996 0.05880 0.05880
1997 0.04620 0.04620
1998 0.04710 0.04710
1999 0.06480 0.06480
2000 0.05720 0.05720
2001 0.07010 0.07010
2002 0.01380 0.01380
2003 0.02030 0.02030
2004 0.01560 0.01560
2005 0.01770 0.01770
2006 0.01770 0.01770
2007 0.01770 0.01770
2008 0.01700 0.01700
2009 0.01840 0.01840

Esd: Additional inter-annual variability in fishing mortality rate. For each historical simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. If this parameter has a positive (non-zero) value, the yearly fishing mortality rate is drawn from a log-normal distribution with a standard deviation (in log space) specified by the value of Esd drawn for that simulation. This parameter applies only to historical projections.

Specified Value(s): 0, 0

Since we have an trend in F this is set to 0.

Historical Effort

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in historical fishing mortality (Esd), with vertical colored lines indicating 3 randomly drawn values used in the time-series plot:

Time-Series

Time-series plot showing 3 trends in historical fishing mortality (OM@EffUpper and OM@EffLower or OM@cpars$Find):

Annual increase in catchability, interannual variability in catchability: qinc, qcv

qinc: Mean temporal trend in catchability (also though of as the efficiency of fishing gear) parameter, expressed as a percentage change in catchability (q) per year. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. Positive numbers indicate an increase and negative numbers indicate a decrease. q then changes by this amount for in each year of the simulation This parameter applies only to forward projections.

Specified Value(s): -0.1, 0.1

The fishery report and assessment provide no compelling reason to expect fishing to become more or less efficient and we set the % annual increase to be very close to zero, -0.1 to 0.1.

qcv: Inter-annual variability in catchability expressed as a coefficient of variation. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This parameter applies only to forward projections.

Specified Value(s): 0.01, 0.03

Given an assessment, variability in fishing efficiency among seasons can be quantified by comparing observed unstandardized catch rates with an index of abundance (or assessed biomass). Here we assume very little variability following a similar analysis for Canary Rockfish: a log-normal standard deviation between 0.01 and 0.03.

Future Catchability

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in fishing efficiency (qcv) and average annual change in fishing efficiency (qinc), with vertical colored lines indicating 3 randomly drawn values used in the time-series plot:

Time-Series

Time-series plot showing 3 trends in future fishing efficiency (catchability):

Fishery gear length selectivity: L5, LFS, Vmaxlen, isRel

L5: Shortest length at which 5% of the population is vulnerable to selection by the gear used in this fleet. Values can either be specified as lengths (in the same units used for the maturity and growth parameters in the stock object) or as a percentage of the size of maturity (see the parameter isRel for more information). For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This value is the same in all years unless cpars is used to provide time-varying selection.

Specified Value(s): 31.44, 38.82

Here we assume that 5% selectivity occurs at age 10. The inverse growth model for the most extreme growth parameters low-K low-Linf / high-K high-Linf predict the specified length range.

LFS: Shortest length at which 100% of the population is vulnerable to selection by the gear used by this fleet. Values can either be specified as lengths (in the same units used for the maturity and growth parameters in the stock object) or as a percentage of the size of maturity (see the parameter isRel for more information). For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This value is the same in all years unless cpars is used to provide time-varying selection.

Specified Value(s): 41.51, 50.11

Uses the same principal above for full selectivity at age 25.

Vmaxlen: Proportion of fish selected by the gear at the asymptotic length (‘’). Upper and Lower bounds between 0 and 1. A value of 1 indicates that 100% of fish are selected at the asymptotic length, and the selection curve is logistic. If Vmaxlen is less than 1 the selection curve is dome shaped. For example, if Vmaxlen is 0.4, then only 40% of fish are vulnerable to the fishing gear at the asymptotic length.

Specified Value(s): 1, 1

Flat-topped selectivity is assumed (as in Canary Rockfish)

isRel: Specify whether selection and retention parameters use absolute lengths or relative to the size of maturity. Single logical value (TRUE or FALSE).

Specified Value(s): FALSE

In this case we are not specifying L5 and LFS as a fraction of length at maturity but rather in absolute units (cm) the same as those of the growth and maturity parameters.

Fishery length retention: LR5, LFR, Rmaxlen, DR

LR5: Shortest length at which 5% of the population is vulnerable to retention by the fleet. Values can either be specified as lengths (in the same units used for the maturity and growth parameters in the stock object) or as a percentage of the size of maturity (see the parameter isRel for more information). For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This value is the same in all years unless cpars is used to provide time-varying selection.

Specified Value(s): 0, 0

Retention follows selectivity.

LFR: Shortest length where 100% of the population is vulnerable to retention by the fleet. Values can either be specified as lengths (in the same units used for the maturity and growth parameters in the stock object) or as a percentage of the size of maturity (see the parameter isRel for more information). For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This value is the same in all years unless cpars is used to provide time-varying selection.

Specified Value(s): 0, 0

Retention follows selectivity.

Rmaxlen: Proportion of fish retained at the asymptotic length (‘’). Upper and Lower bounds between 0 and 1. A value of 1 indicates that 100% of fish are retained at the asymptotic length, and the selection curve is logistic. If Rmaxlen is less than 1 the retention curve is dome shaped. For example, if Rmaxlen is 0.4, then only 40% of fish at the asymptotic length are retained.

Specified Value(s): 1, 1

Retention follows selectivity.

DR: Discard rate, defined as the proportion of fully selected fish that are discarded by the fleet. Upper and Lower bounds between 0 and 1, with a value of 1 indicates that 100% of selected fish are discarded. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided.

Specified Value(s): 0, 0.1

The general discard rate. What fraction of fish across all size and age classes are discarded? There is general discarding across all size-classes in some fisheries. We assume that this is up to 10% in the Canary Rockfish fishery.

Current Year: CurrentYr

CurrentYr: The last historical year simulated before projections begin. Single value. Note that this should match the last historical year specified in the Data object, which is usually the last historical year for which data is available.

Specified Value(s): 2009

Existing Spatial Closures: MPA

MPA: Logical argument (TRUE or FALSE). Creates an MPA in Area 1 for all years if true is selected. Defaults to FALSE.

Specified Value(s): FALSE

No justification provided.

Obs Parameters

The observation model parameter are taken from the Generic_Obs model subject to a few addtional changes which are documented here.

Catch statistics: Cobs, Cbiascv, CAA_nsamp, CAA_ESS, CAL_nsamp, CAL_ESS

Cobs: Observation error around the total catch. Observation error in the total catch is expressed as a coefficient of variation (CV). Cobs requires upper and lower bounds of a uniform distribution, and for each simulation a CV is sampled from this distribution. Each CV is used to specify a log-normal error distribution with a mean of 1 and a standard deviation equal to the sampled CV. The yearly observation error values for the catch data are then drawn from this distribution. For each time step the simulation model records the true catch, but the observed catch is generated by applying this yearly error term (plus any bias, if specified) to the true catch.

Specified Value(s): 0.1, 0.2

Catches are observed more precisely than the Generic_Obs object with a CV of between 10 and 20 per cent.

Cbiascv: Log-normally distributed coefficient of variation controlling the sampling bias in observed catch for each simulation. Bias occurs when catches are systematically skewed away from the true catch level (for example, due to underreporting of catch or undetected illegal catches). Cbiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years.

Specified Value(s): 0.05

Mean bias (under / over reporting) in catches is assumed to be small with a CV of 0.05 95% of simulations are reported between 90% and 110% of true simulated catches.

CAA_nsamp: Number of catch-at-age observations collected per time step. For each time step a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. Positive integers.

Specified Value(s): 100, 200

An initial guess at 200 per year

CAA_ESS: Effective sample size of catch-at-age observations collected per time step. For each time step a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. CAA_ESS should not exceed CAA_nsamp. Positive integers.

Specified Value(s): 25, 50

Effective sample size is identical to CAA_nsamp

CAL_nsamp: Number of catch-at-length observations collected per time step. For each time step a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. Positive integers.

Specified Value(s): 100, 200

Assumed to be the same as catch at age

CAL_ESS: Effective sample size. For each time step a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. CAL_ESS should not exceed CAL_nsamp. Positive integers.

Specified Value(s): 25, 50

Assumed to be the same as catch at age

Index imprecision, bias and hyperstability: Iobs, Btobs, Btbiascv, beta

Iobs: Observation error in the relative abundance index expressed as a coefficient of variation (CV). Iobs requires upper and lower bounds of a uniform distribution, and for each simulation a CV is sampled from this distribution. Each CV is used to specify a log-normal error distribution with a mean of 1 and a standard deviation equal to the sampled CV. The yearly observation error values for the index of abundance data are then drawn from this distribution. For each time step the simulation model records the true change in abundance, but the observed index is generated by applying this yearly error term (plus any bias, if specified) to the true relative change in abundance. Positive real numbers.

Specified Value(s): 0.1, 0.25

Relative abundance indices are assumed to be observed imprecisely CV of between 20 and 30 per cent.

Btobs: Observation error in the absolute abundance expressed as a coefficient of variation (CV). Btobs requires upper and lower bounds of a uniform distribution, and for each simulation a CV is sampled from this distribution. Each CV is used to specify a log-normal error distribution with a mean of 1 and a standard deviation equal to the sampled CV. The yearly observation error values for the absolute abundance data are then drawn from this distribution. For each time step the simulation model records the true abundance, but the observed abundance is generated by applying this yearly error term (plus any bias, if specified) to the true abundance. Positive real numbers.

Specified Value(s): 0.2, 0.5

As Generic_Obs.

Btbiascv: Log-normally distributed coefficient (CV) controlling error in observations of the current stock biomass. Bias occurs when the observed index of abundance is is systematically higher or lower than the true relative abundance. Btbiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0.33

The bias in the absolute abundance index is assumed to be reasonably high and could be 1/3 to 3 times the true value.

beta: A parameter controlling hyperstability/hyperdepletion in the measurement of abundance. For each simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. Values below 1 lead to hyperstability (the observed index decreases more slowly than the true abundance) and values above 1 lead to hyperdepletion (the observed index decreases more rapidly than true abundance). Positive real numbers.

Specified Value(s): 0.66, 1.5

Since the survey is carried out according to a systematic design we assume that it varies roughly proportionally to real abundance and specify a beta parameter between 2/3 and 3/2.

Bias in maturity, natural mortality rate and growth parameters: LenMbiascv, Mbiascv, Kbiascv,t0biascv, Linfbiascv

LenMbiascv: Log-normal coefficient of variation for sampling bias in observed length at 50 percent maturity. LenMbiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0.1

As Generic_Obs.

Mbiascv: Log-normal coefficient of variation for sampling bias in observed natural mortality rate. Uniform distribution lower and upper bounds. Mbiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0.2

As Generic_Obs.

Kbiascv: Log-normal coefficient of variation for sampling bias in observed growth parameter K. Kbiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0.05

Twice as accurate as Generic_Obs with K values within plus or minus 5% of true value.

t0biascv: Log-normal coefficient of variation for sampling bias in observed t0. t0biascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0

We choose not to simulate bias in this growth parameter and assume in all cases it is correct.

Linfbiascv: Log-normal coefficient of variation for sampling bias in observed maximum length. Linfbiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0.03

Twice as accurate as Generic_Obs with L-infinity values within plus or minus 5% of true value.

Bias in length at first capture, length at full selection: LFCbiascv, LFSbiascv

LFCbiascv: Log-normal coefficient of variation for sampling bias in observed length at first capture. LFCbiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0.05

Given the reasonably extensive length sampling data, it is straightforward to estimate Length at First Capture for rockfish from the length frequency data and this is likely to be reasonably well known without substantial bias.

LFSbiascv: Log-normal coefficient of variation for sampling bias in length-at-full selection. LFSbiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0.05

As Generic_Obs

Bias in fishery reference points, unfished biomass, FMSY, FMSY/M ratio, biomass at MSY relative to unfished: FMSY_Mbiascv, BMSY_B0biascv

FMSY_Mbiascv: Log-normal coefficient of variation for sampling bias in estimates of the ratio of the fishing mortality rate that gives the maximum sustainable yield relative to the assumed instantaneous natural mortality rate. FMSY/M. FMSY_Mbiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0.15

A number of MPs aim to fish at a fixed rate proportional to the estimate of M (e.g. Fratio). Other MPs use this ratio to undertake stock reduction analysis (e.g. DB-SRA). Given the references above we set this to be moderately inaccurate given a CV of 0.15.

BMSY_B0biascv: Log-normal coefficient of variation for sampling bias in estimates of the BMSY relative to unfished biomass (BMSY/B0). BMSY_B0biascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0.05

We have assigned a relatively precise CV for potential accuracy at 0.05.

Management targets in terms of the index (i.e., model free), the total annual catches and absolute biomass levels: Irefbiascv, Crefbiascv, Brefbiascv

Irefbiascv: Log-normal coefficient of variation for sampling bias in the observed relative index of abundance (Iref). Irefbiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0.2

Here we assume that the index and MSY (a desirable catch level) can be known more accurately than a desirable absolute biomass level (e.g. BMSY) and assign these a range determined by a CV of 0.2.

Crefbiascv: Log-normal coefficient of variation for sampling bias in the observed reference catch (Cref). Crefbiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0.2

As Irefbiascv.

Brefbiascv: Log-normal coefficient of variation for sampling bias in the observed reference biomass (Bref). Brefbiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0.5

Arbitrarily we make this twice as potentially biased as Iref and Cref.

Depletion bias and imprecision: Dbiascv, Dobs

Dbiascv: Log-normal coefficient of variation for sampling bias in the observed depletion level. Dbiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0.5

Here we assign an arbitrary value of 0.25 which is relatively imprecise and means that assume depletion could up to double or half of the true simulated value.

Dobs: Log-normal coefficient of variation controlling error in observations of stock depletion among years. Observation error in the depletion expressed as a coefficient of variation (CV). Dobs requires the upper and lower bounds of a uniform distribution, and for each simulation a CV is sampled from this distribution. Each CV is used to specify a log-normal error distribution with a mean of 1 and a standard deviation equal to the sampled CV. The yearly observation error values for the depletion data are then drawn from this distribution. For each time step the simulation model records the true depletion, but the observed depletion is generated by applying this yearly error term (plus any bias, if specified) to the true depletion.

Specified Value(s): 0.05, 0.1

In a data-limited situation it is unlikely that radically new data would become available regarding depletion meaning that while estimates may be biased, they are likely to be relatively precise. We assign a level of imprecision consistent with observations of catch rate data among years at between 0.05- 0.1.

Recruitment compensation and trend: hbiascv, Recbiascv, sigmaRbiascv

hbiascv: Log-normal coefficient of variation for sampling persistent bias in steepness. hbiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0.2

The stock assessment provides very little support for particular values of recruitment compensation. In DLMtool this is parameterized as steepness (the fraction of unfished recruitment at 20% of unfished spawning biomass, a value ranging from 0.2-1). Here we assume that any MP could get this wrong by a large margin.

Recbiascv: Log-normal coefficient of variation for sampling persistent bias in recent recruitment strength. Recbiascv requires the upper and lower bounds of a uniform distribution, and for each simulation a CV is sampled from this distribution. Each CV is used to specify a log-normal error distribution with a mean of 1 and a standard deviation equal to the sampled CV. The yearly bias values for the depletion data are then drawn from this distribution. Positive real numbers.

Specified Value(s): 0.1, 0.2

As Generic_Obs

sigmaRbiascv: Log-normal coefficient of variation for sampling persistent bias in recruitment variability. sigmaRbiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years. Positive real numbers.

Specified Value(s): 0.1

No justification provided.

Effort: Eobs, Ebiascv

Eobs: Observation error around the total effort. Observation error in the total effort is expressed as a coefficient of variation (CV). Eobs requires upper and lower bounds of a uniform distribution, and for each simulation a CV is sampled from this distribution. Each CV is used to specify a log-normal error distribution with a mean of 1 and a standard deviation equal to the sampled CV. The yearly observation error values for the effort data are then drawn from this distribution. For each time step the simulation model records the true effort, but the observed effort is generated by applying this yearly error term (plus any bias, if specified) to the true effort.

Specified Value(s): 0, 0

No justification provided.

Ebiascv: Log-normally distributed coefficient of variation controlling the sampling bias in observed effort for each simulation. Bias occurs when effort is systematically skewed away from the true effort level. Ebiascv is a single value specifying the standard deviation of a log-normal distribution with a mean of 1 and a standard deviation equal to the sampled CV. For each simulation a bias value is drawn from this distribution, and that bias is applied across all years.

Specified Value(s): 0

No justification provided.

Obs Plots

Observation Parameters

Catch Observations

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in catch observations (Csd) and persistent bias in observed catch (Cbias), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Time-series plots of catch observation error for historical and projection years:

Depletion Observations

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in depletion observations (Dobs) and persistent bias in observed depletion (Dbias), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Time-series plots of depletion observation error for historical and projection years:

Abundance Observations

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in abundance observations (Btobs) and persistent bias in observed abundance (Btbias), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Time-series plots of abundance observation error for historical and projection years:

Index Observations

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in index observations (Iobs) and hyper-stability/depletion in observed index (beta), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Time-series plot of 3 samples of index observation error:

Plot showing an example true abundance index (blue) with 3 samples of index observation error and the hyper-stability/depletion parameter (beta):

Recruitment Observations

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in index observations (Recsd) , with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Timeseries plots of observeration error for recruitment:

Composition Observations

Sampled Parameters

Histograms of 48 simulations of catch-at-age effective sample size (CAA_ESS) and sample size (CAA_nsamp) and catch-at-length effective (CAL_ESS) and actual sample size (CAL_nsamp) with vertical colored lines indicating 3 randomly drawn values:

Parameter Observations

Sampled Parameters

Histograms of 48 simulations of bias in observed natural mortality (Mbias), von Bertalanffy growth function parameters (Linfbias, Kbias, and t0bias), length-at-maturity (lenMbias), and bias in observed length at first capture (LFCbias) and first length at full capture (LFSbias) with vertical colored lines indicating 3 randomly drawn values:

Reference Point Observations

Sampled Parameters

Histograms of 48 simulations of bias in observed FMSY/M (FMSY_Mbias), BMSY/B0 (BMSY_B0bias), reference index (Irefbias), reference abundance (Brefbias) and reference catch (Crefbias), with vertical colored lines indicating 3 randomly drawn values:

Imp Parameters

Output Control Implementation Error: TACFrac, TACSD

TACFrac: Mean fraction of recommended TAC that is actually taken. For each historical simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This value is the mean TAC fraction obtained across all years of that simulation, and a yearly TAC frac is drawn from a log-normal distribution with the simulation mean and a coefficient of variation specified by the value of TACSD drawn for that simulation. If the value drawn is greater than 1 the amount of catch taken is greater than that recommended by the TAC, and if it is less than 1 the amount of catch taken is less than that recommended by the TAC. Positive real numbers.

Specified Value(s): 1, 1.1

Here we assume that the actual catches can be up to 10% higher than the recommended TAC.

TACSD: Log-normal coefficient of variation in the fraction of recommended TAC that is actually taken. For each historical simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This value is used, along with the TACFrac drawn for that simulation, to create a log-normal distribution that yearly values specifying the actual amount of catch taken are drawn from. Positive real numbers.

Specified Value(s): 0.05, 0.1

We assume that the bias in the actual catch is relatively consistent between years and set the range for this parameter to a low value.

Effort Control Implementation Error: TAEFrac, TAESD

TAEFrac: Mean fraction of recommended TAE that is actually taken. For each historical simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This value is the mean TAE fraction obtained across all years of that simulation, and a yearly TAE frac is drawn from a log-normal distribution with the simulation mean and a coefficient of variation specified by the value of TAESD drawn for that simulation. If the value drawn is greater than 1 the amount of effort employed is greater than that recommended by the TAE, and if it is less than 1 the amount of effort employed is less than that recommended by the TAE. Positive real numbers.

Specified Value(s): 1, 1.1

We have little information to inform this parameter, and set the implementation error in effort equal to the TAC implementation error.

TAESD: Log-normal coefficient of variation in the fraction of recommended TAE that is actually taken. For each historical simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This value is used, along with the TAEFrac drawn for that simulation, to create a log-normal distribution that yearly values speciying the actual amount of efort employed are drawn from. Positive real numbers.

Specified Value(s): 0.05, 0.1

We assume that the bias in the effort is relatively consistent between years and set the range for this parameter to a low value.

Size Limit Control Implementation Error: SizeLimFrac, SizeLimSD

SizeLimFrac: Mean fraction of recommended size limit that is actually retained. For each historical simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This value is the mean size limit fraction obtained across all years of that simulation, and a yearly size limit fraction is drawn from a log-normal distribution with the simulation mean and a coefficient of variation specified by the value of SizeLimSD drawn for that simulation. If the value drawn is greater than 1 the size of fish retained is greater than that recommended by the size limit, and if it is less than 1 the amount of size of fish retained is less than that recommended by the size limit. Positive real numbers.

Specified Value(s): 1, 1.1

We assume that, on average, a size limit would be well-implemented.

SizeLimSD: Log-normal coefficient of variation in the fraction of recommended size limit that is actually retained. For each historical simulation a single value is drawn from a uniform distribution specified by the upper and lower bounds provided. This value is used, along with the SizeLimFrac drawn for that simulation, to create a log-normal distribution that yearly values speciying the actual fraction of the size limit retained are drawn from. Positive real numbers.

Specified Value(s): 0.05, 0.1

We assume that the implementation of the size limit is relatively consistent between years.

Imp Plots

Implementation Parameters

TAC Implementation

Sampled Parameters

Histograms of 0 simulations of inter-annual variability in TAC implementation error (TACSD) and persistent bias in TAC implementation (TACFrac), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Time-series plots of 0 samples of TAC implementation error by year:

TAE Implementation

Sampled Parameters

Histograms of 0 simulations of inter-annual variability in TAE implementation error (TAESD) and persistent bias in TAC implementation (TAEFrac), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Time-series plots of 0 samples of TAE implementation error by year:

Size Limit Implementation

Sampled Parameters

Histograms of 0 simulations of inter-annual variability in size limit implementation error (SizeLimSD) and persistent bias in size limit implementation (SizeLimFrac), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Time-series plots of 0 samples of Size Limit implementation error by year:

Historical Simulation Plots

Historical Time-Series

Spawning Biomass

Depletion

Time-series plots of SB/SB0:

Absolute

Time-series plots of absolute SB:

Vulnerable Biomass

Depletion

Time-series plots of VB/VB0:

Absolute

Time-series plots of absolute VB:

Total Biomass

Depletion

Time-series plots of B/B0:

Absolute

Time-series plots of absolute B:

Recruitment

Relative

Time-series plot of recruitment relative to R0:

Absolute

Time-series plot of absolute recruitment:

Catch

Relative

Time-series of catch relative to the current year:

Absolute

Time-series of absolute catch:

Historical Fishing Mortality

Historical Time-Series

Time-series of historical fishing mortality: